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Q ; Abstract 



The stability of a Fermi ball (F-ball), which is a kind of non-top ological soliton accompanying 
the breakdown of the approximate Z2 symmetry, is investigated in three situations: the case it is 



Qh, electrically neutral, the case it is electrically charged and unscreened, and the case it is electrically 



I charged and screened. We argue only the third case is physically meaningful since the neutral 



F-ball is unstable and the case of the unscreened charged one is observationally excluded when 
^ . it has a sizable contribution to CDM. We found that the energy scale of the breakdown of the 

. . , approximate Z2 symmetry v should satisfy v < 3 x 10 GeV if the F-ball is a main component of 

CDM. 

PACS numbers: 05.45.Yv, 95.35. +d 
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I. INTRODUCTION 



An F-ball, a kind of non-topological soliton, was introduced as a candidate for CDM 
|l|, |2[. A similar object is also suggested to play an important role in baryogenesis 0. 

The simplest type of the F-ball is a bubble of false vacuum surrounded by a thin domain 
wall, on which many zero-mode fermions are attached. The surface tension of the wall is 
balanced with the Fermi pressure due to the zero-mode fermions in the wall. The F-ball 
may be stabilized owing to this balance between the surface energy and the Fermi energy. 
Though the F-ball is stable against the variation of the radius in case it keeps a spherical 
shape, it may be unstable against deformation from sphere. Macpherson and Campbell 
pointed out in their pioneering work [|l| that such F-balls are unstable against deformation 
from the spherical shape, and that they should finally fragment into a number of tiny F-balls 
|T8| (see Fig.|l]for illustration of a fragmenting F-ball). Suppose there are a large number of 



J 



(a) (b) (c) (d) 

FIG. 1: A deformation and a fragmentation of the F-ball. (a) Sphere, (b) Cigar-like deformation, 
(c) Shape narrow in the middle, (d) Fragmented F-balls. 

such tiny F-balls in the present universe contributing sizably to CDM. The cosmic flux of 
these objects should be large, proportional to the inverse of the mass, and they should have 
been observed unless the scattering cross section with ordinary matter is very small. Here, 
we are not interested in such a new kind of weakly interacting particles since we need too 
many new assumptions in introducing new tiny particles and have only indirect means to 
detect such particles M. 



Morris then introduced a new idea, an electrically charged F-ball which is considered 
to be stabilized due to the repulsive long-range Coulomb force . Since the charged F-bahs 
can be large and heavy in this case, they can sizably contribute to CDM without contradict- 



ing the present observations of the cosmic flux. These F-balls with a large cross section of 
interaction with ordinary matter are interesting because of their future detectabihty. (We 
nevertheless note that the cross section per unit mass for the F-ball a /Mf is not so large as 
the one which was recently proposed ^ for the self-interacting dark matter. We do not 
go into details since the above issue is still controversial (see the discussions in Sec.|V|).) 

In the present paper, we flrst examine the instability of the neutral F-ball against the 
deformation from the spherical shape. It is known [|I|] that the spherical shape is unstable 
against deformation due to the flnite volume energy within the thin-wall approximation. We 
investigate the next-to-leading order approximation, the curvature effect of the wall, since it 
has a shape- dependent contribution to the total energy, and is important when the volume 
energy is very small. We flnd that the neutral F-ball is unstable even in the presence of the 
curvature effect. 

We next examine whether the F-ball is stable when it is electrically charged and un- 
screened. Since the total energy of such an F-ball is higher than the sum of the energies of 
the fragmented F-balls, it fragments into smaller F-balls in finite time. Its lifetime is, how- 
ever, longer than the age of the universe if there is an energy barrier high enough between 
the state of an F-ball and that of fragmented F-balls. We see that the Coulomb force does 
make the lifetime of the F-ball long enough. 

We thirdly examine whether the F-ball has a lifetime long enough to survive until present 
if it is screened due to the electrons or positrons in the thermal bath of the early universe. 
Though the long ranged Coulomb force becomes short ranged in the thermal bath, the 
lifetime of the F-ball is found long enough in the proper region of the parameters. 

We finally examine which kind of F-balls satisfy the conditions for a main component of 
CDM. We find that only the screened F-ball meets the conditions after taking into account 
the stability of the F-ball, observational constraints, and cosmological considerations. 

The above contents are organized as follows: We first consider electrically neutral F-balls 
and discuss their stability in Sec.|I[ Small F-balls which are electrically charged and not 



screened by electrons are considered in Sec. III. Screened ones are considered in Sec.IV 



this is the main theme in the present paper. In Sec.|V|, we discuss certain constraints on 
the parameters obtained from various conditions. We have summaries and discussions in 
SeclY^. The detailed calculations are given in the Appendices. 
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II. STABILITY OF ELECTRICALLY NEUTRAL F-BALL 



We briefly introduce the electrically neutral F-ball proposed by Macpherson and Campbell 
pi] to make clear the notations and technical terms. The Lagrangian density is given by 



(1) 



where and ipj are a scalar and a fermion field, respectively, and G is a Yukawa coupling 
constant. Here, f/(0) is an approximate double-well potential [ pO[] , 



(2) 



The first term has the Z2 symmetry under ^ —(f), and the second term violates the 
symmetry, though we assume it is much smaller than the first one (see Fig.|3). Thus, the 
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FIG. 2: The biased potential for A = 1. Two almost degenerate vacua witli the energy difference 
of e exist: tlie true vacuum with 6 = v and the false vacuum with 6 = —v. 



Lagrangian density has an approximate Z2 symmetry, which we call 'biased Z2 symmetry' 
[|1|] hereafter. 

After the phase transition breaks spontaneously the biased Z2 symmetry at T = Tph ~ v, 
there arise two almost degenerate vacua with the energy density difference of e. They are 
the true vacuum with cj) = v and the false one with = —v. Then a domain wall is 
produced between two vacua, and the fermions are captured as zero-modes in the domain 
wall Though the effect of e should be negligible soon after the phase transition, it 

gets more important as the universe expands. If the volume energy becomes larger than the 
surface energy during the expansion, the true vacuum is energetically favored and enlarges 
its volume. At the same time, the region of false vacuum gets diminished, and becomes a 



small confined region surrounded by the domain wall. Such a region of false vacuum would 
continue to shrink due to the surface tension and the volume effect proportional to e, and 
finally disappear if we neglect the Fermi pressure caused by the zero-mode fermions trapped 
in the wall. In our case, however, the Fermi pressure is not negligible and stops the shrinkage 
of the region, when it gets balanced with the surface tension. Such a bubble of the false 
vacuum with zero-mode fermions trapped in the surrounding wall is called F-ball We 
take the temperature for the F-ball production as T = Tf. Because Ty is somewhat lower 
than Tph, Tf is assumed to be Ty < O.lf in the present paper. 

We consider a thin walled F-ball, the size of which is much larger than the thickness of 
the domain wall. For simplicity, we neglect the small energy difference e and take the exact 
Z2 symmetry until Sec.0. Then the Lagrangian density is 

^ = lid,<PY + ^Pfii^^d, - G<P)i,f - ^ - ■ (3) 

The total energy of an F-ball is expressed as 

Etot = Es + Ef , (4) 

where Eg and Ej are a contribution from the domain wall and that from the zero-mode 
fermions, respectively. In the thin wall case, Fermi gas in the wall distributes two dimen- 
sionally on the surface of the F-ball p2[. Taking w-axis normal to the surface which is 
positive (negative) for the outside (inside) of the F-ball, we express the fermion number 
density as 

Uf = af6{w) , (5) 
with aj the surface density of the fermion number, and the fermi energy, 

= 1^ / dS af . (6) 

We obtain the surface energy. 



where Ri and R2 are the radii of principal curvature of the F-ball surface (see Fig.^ and 
Appendix 0). 

If the fermion-scalar coupling G is not too small, the F-ball is stable against releasing a 
fermion with mass ~ Gv from the wall to the vacuum. The number of the fermions of the 



5 



FIG. 3: The local orthogonal coordinates {u,v,w), the origin of which is on the surface of the 
F-ball. We take u and v along the lines of curvature on the surface, and w normal to the surface. 
Here, Ri {R2) is a radius of principal curvature with respect to u = (u = 0) at the origin. 




(8) 



is then conserved. Thus, the energy of the F-ball is obtained by minimizing 




(9) 



with n the Lagrange multiplier. We first minimize it with respect to o"/. Noting that Eg is 
independent of af, we obtain the uniform distribution of the fermions, 

^f^i-^ (10) 

with S the area of the F-ball surface. We next minimize the energy with respect to 0. 
Assuming that the variation of in the direction parallel to the F-ball surface is much 
smaller than that in the direction w, we get the Euler-Lagrange equation for 0, 




(11) 



with the symbol 'dot' in or standing for the derivative with respect to w. Here, fulfills 
the boundary condition. 



= 



+V {w > Sn) 

= 0) 

-V {w < -5n) 



(12) 
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with Sn the width of the wall. We take the leading contribution, (po, in the thin- wall expansion 
with Ri, R2 ^ oo. It satisfies 



A 



'o(0O 



(13) 



with 



+V {W > 5n) 

(w = 0) 



(14) 



The solution for 



^0, 



V tanh 



w 

5n 



with 6n = 2/ (\/Af ), minimizes the energy under the condition, Eq.(|ToD, to be 



(0) 
tot 



S5 + 



(15) 



(16) 



with Tj = 2^/Xv^/3 the surface tension in the wall. We finally minimize E^^j with respect to 
S, and obtain 

E^'J = {12nJ:y/''Nf , (17) 

with the of the surface, S = {2^ /"SLf/^Nf. An F-ball can deform with this surface 
area being kept constant (see Fig.|I](a) to Fig.|l|(c)). Moreover, it has the same energy even 
if it splits into some smaller F-balls (see Fig.|l|(d)) since Ef^^ is proportional to the number 
of fermions. Thus, we cannot tell whether the F-ball is stable or not against deformation 
and fragmentation into pieces, within the thin-wall approximation which corresponds to 
the leading order in the thin-wall expansion. This is the same result as what Macpherson 
and Campbell derived in their paper They concluded further that the neutral F-ball 
is unstable, taking into account the volume effect proportional to e. We here consider the 
next-to-leading order contribution in the thin-wall expansion, i.e., the curvature effect, since 
such effect becomes important when it is comparable to the volume energy effect. 
We keep the terms up to the next-to-leading order for 0, 



+ 01 



(18) 



with 01 smaller than 0o by the order of 5„/-R 

A 
2 



From Eqs. ([TI|) and ([12|), 0i satisfies 



1 ^ 1 



^0 , 



(19) 



(20) 



and the boundary condition, 

{\W\ > 6n) 

{w = 0) . 

The analytic solution ||2^ for (pi is given in Appendix (see Fig.^] for the illustration of (po 
and (pi). Substituting the solution for (pi and that for (pQ into Etot and keeping the terms up 
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FIG. 4: The perturbative solutions for the domain wall field (p. Here, 0o is the leading term and 
(pi the next-to- leading term in the thin- wall expansion. Note that (pQ is rescaled as (po/v with v 
the symmetry breaking scale, and (pi as VX(pi with (pi = 2RiR2(pi/{Ri + R2) {Ri and R2 are the 
radii of principal curvature of the F-ball surface). 



to the order of 0{{6J R)^Ei^l) we get (see Appendix @) 



p _ p(0) , p(2) 



(21) 



where 



eI^I = Cn I dS + Dr, I dS 



\R1R2 



1 1 \2 \R1R2 



Ri R2J RiR 



Here, C„ and Dn are given by 



■ 2 



dw i^W 00 - 0001 

0.25v 



(22) 



(23) 



with 01 = 2RiR2(j)i/{Ri + R2). 

First, we examine the stabihty of the spherical F-ball against the perturbative deformation 
from the spherical shape (see Fig.|(b)). For this purpose, we take the curvature radius on 
the F-ball surface positive. In this case, the first term in Eq.(^2]) is constant according to 



the Gauss-Bonnet theorem |£6|. Thus, the shape-dependent sector in Etot comes from the 
second term in Eg . (^2]) . Because it is semi-negative definite, and zero only for the spherical 
shape, the spherical F-ball is unstable against the perturbative deformation. 

We next consider the stability of the F-ball against fragmentation into smaller ones (see 
Fig.|l](d)). In this case, the fragmentation does not change i^^i because it only depends on 
the fermion number. Here, we compare the energy of the spherical F-ball with that of the 
fragmented spherical F-balls for simplicity. We see that the second term in Eq.(p^) vanishes. 
Thus, only the first term in Eg. (|22|) changes in the fragmentation. Since it is proportional 
to the number of F-balls due to the Gauss-Bonnet theorem, its negative sign shows that 
the fragmentation is energetically promoted. Therefore, we find that the F-ball is unstable 
against not only the deformation from the spherical shape but also the fragmentation into 
smaller ones. 

As pointed out by Macpherson and Campbell, the fragmentation will continue until the 
thin-wall expansion gets invalid. If the fragmented tiny F-balls sizably contribute to CDM, 
there should be a large number of them and a large cosmic flux of them as well. The present 
observations for the dark matter search tell that they can interact with ordinary matter only 
very weakly. Such F-balls are not considered here since their detectability was discussed in 
Ref.e. 



III. STABILITY OF ELECTRICALLY CHARGED F-BALL IN THE UN- 
SCREENED CASE 

In the present section, we consider the electrically charged F-ball proposed by Morris , 
in which the fermions trapped on the surface have electric charge. Assigning the electric 
charge of +e to the fermion for simplicity, we consider an F-ball carrying the electric charge 
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+eNf. The energy of the F-ball is expressed as 

Etot = E, + Ef + E,, (25) 

with Ec the Coulomb potential energy. Since the Coulomb energy is proportional to the 
fermion number squared Ec oc Nj, the total energy per unit fermion number Etot/Nf is 
an increasing function of Nf. Thus, when the F-ball fragments into smaller ones, the total 
energy apparently decreases. We therefore cannot conclude that the electric force stabilizes 
the F-ball from fragmentation only by comparing the energies of the two states. Morris, 
however, suggested that the inclusion of the long ranged repulsive force helps the F-ball 
to be stabilized. Since it is not so clear whether the electric repulsive force stabilizes the 
F-ball or not, we here verify it using the perturbative method. 

For this purpose, we investigate whether the energy barrier exists between the F-ball 
before the fragmentation and that after it. Let us consider a spherical F-ball and its per- 
turbative deformation from the spherical shape. Assuming the rotational invariance around 
z-axis for simplicity, we express the surface of the F-ball in the polar coordinates as 

^j{e,^)=r{9)e, (26) 



with e a unit vector. 



We expand r{9) as 



e = (sin 6 cos (f, sin 6 sin (f, cos 6) . (27) 



rie) = R\^l + J2aiPi{cose)j 

= R{1 + Sr) , (28) 

where a^'s are small coefficients of the perturbative expansion, and P/(cos^) is the Legendre 
polynomial. We write the surface element as 

dS = R^{l + 6S{ai))dQ , (29) 

with dQ = sin 9d9dip. We also expand the surface density of the fermion as 

^ 'l + 6af) , (30) 



47rP2(i + SS) 
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where q's are small expansion coefficients. Note that the above expression satisfies the 
condition for the fermion number Nf to be conserved (see Eq.(|^)). 

Keeping the terms up to the second power of a; and q (see Appendix 0), we estimate, 
within the thin-wall approximation, the surface energy. 



4vrE/?Ml + ET7r^ , (31) 



"2/ + l^ 



the fermi energy. 



3/2 



and the Coulomb potential energy, 

SnJ J 



where 



(l + ^) «f , (34) 

Hf ^ -i (l . '1±11) al . hi , (35) 

+ 2 2(/-l) , 1 ,2 

h^P = ; of -afii + — bj , (36) 

2+1 ' 2+12+1'' ^ ^ 



with bi = ci — 2a;. Using the effective radius. 



we express the total energy as 



where 



Ei,i = 47rSi?2 + ^l!^ 1 + V^^ +^ 1 + V^^^ , (38) 



= , (39) 

<7« = (/ - 1) (- , , + - - + • (40) 

^ ' \ 2(2/ + 1) 4; ' 2/ + 1 ' ' 2/ + 1 ^ ^ 
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Minimizing Etot with respect to Re, we have 

fo / \ 2/3 
1 + ^^ 



3(1 + ^^)'=^ 



with the effective radius, 

/Nf 

Re 




Since the last term in the brackets { } in Eq.(^) is semi-positive definite and zero 
only for bi = and ai = bi = for / > 2, Efot takes local minimum when the F-ball is 
spherical (we note that the parameter ai corresponds to the small parallel displacement 
and not to the deformation of the F-ball). Thus, the spherical F-ball is stable against the 
perturbative deformation from the spherical shape. It leads to the existence of the energy 
barrier and should suppress the fragmentation of the spherical F-ball. If the suppression 
of the fragmentation rate is strong, the F-balls would be effectively stable in the universe 
(for the detailed discussion on their lifetime, see Appendix 0) and could be the candidates 
for CDM. However, if Nf is very large, the F-ball should have such a strong electric field 
that should be quantum-mechanically screened by electrons or positrons with the screening 
length much shorter than the F-ball radius, as we see in the next section. Therefore, the 
argument of this section cannot be applied for the F-ball with very large Nf. 

The argument of this section is applicable only for the F-balls with small Nf, which have 
the radius much smaller than the screening length. Such F-balls should give too large cosmic 
flux to be compatible with the present observations if the F-balls contribute sizably to CDM 
(the more details are given in Sec.0 and Ref.|§]). We therefore find that the unscreened 
charged F-ball cannot be a main component of CDM. 
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IV. STABILITY OF ELECTRICALLY CHARGED F-BALL IN THE SCREENED 
CASE 



In this section, we examine the stabihty of such a large F-ball that is electrically charged 
and screened due to the electrons produced quantum mechanically or thermally. We call the 
fermion trapped on the F-ball surface as a 'heavy fermion' in order to distinguish it from 
the electron; the heavy fermion is almost massless in the two dimensional surface while it is 
as heavy as Gv in the true vacuum. The property of the screened F-ball is quite different 
from that of unscreened one, because the long ranged electric force becomes short ranged as 
a result of the screening. Let us investigate the screening effect by introducing the Coulomb 
potential energy Vg for electron (— K for positron) with = —cAq. We use the Thomas- 
Fermi method 10| to deal with this problem [^. The Helmholtz free energy is expressed 



as 

Ftot = Fn + , (43) 

where F„ is the non-Coulombic term that is the same as the energy for the neutral F-ball 
(see Sec.||). The second term Fg^ is the contribution from the screening electrons, 

dSdw (1 + ;^)(1 + -j^^) |-^(VV;)' + - njV^j , (44) 



with jFg the free energy density of electron and positron in the Coulomb potential |[Tl[| , 

= / d^P {log (l + e-^) + log (l + e-^) } . (45) 



Here, the kinetic energy E is given hy E = y + with mg the electron mass, and the 
density of the heavy fermion nj is given by Uf = af6{w). (Note that the negative sign of 
the first term in { } in Eq.(^4D gives the correct electric energy density, -|-E^/2, after the 
free energy is extremized with respect to Vg.) Taking into account the conservation of the 
heavy fermion number, we obtain the free energy by extremizing Ffot + f^{Nf — J dS cxf) 
with respect to 0, ct/ and V^. The procedure of extremization with respect to (p and o"/ is 
the same as in Sec.0, and we refer to the discussions and results there. 

Now let us extremize the free energy with respect to Ve- Assuming that the variation of 
Ve in the direction parallel to the F-ball surface is much smaller than that in the normal 
direction, we get the Thomas-Fermi equation, 

^"^e + ^ f + -^^) Ve = rifiw) - n,iw) , (46) 



\Ri + w R2 + W 
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where the symbol 'dot' stands for the partial derivative with respect to w. Here, ne{w) = 
dJ-'e/dVe is the expectation value of the difference between the electron density and the 
positron density. The second term of l.h.s. in Eq.(^6|) comes from the curvature effect, 
which was not taken into account in the previous works. We here impose the boundary 
condition, 



for \w\ 



(47) 



with 6sc the typical screening length. For the case, T 3> rrie, the electron number density is 
expressed as 



- - a o 2 • (48) 

Note that is positive since we take Vg negative, assuming the F-ball carrying the positive 
electric charge. In the first place, let us obtain the leading contribution in the thin-wall 
expansion by taking -Ri, -R2 00. We get 



sinh (^) 



(49) 



where 



eT ' 

Arcosh"^ 



(50) 



V2nT^ 



2n^T^ + 3e V2 



2V6' 



vr 



> T) . 



(51) 



Substituting V^^^ into Fsc and keeping the terms of the leading order, we get the leading 
order of Fgc as 

1 



dS I dw 



2 ^ ^0 



6 



127r2 



- crfVo6{ 



w 



dS COth^ - -COth + - 



A 



2^2' 



Tie 



\ 3^ 



dS 



a 



3/2 



2 

V>T) 



A 



(52) 



From the above equation, we see that the heavy fermions distribute uniformly over the 
surface, as is the case with Eq.(|TO|). Then, in Eq.(0) and V^^'^ in Eq.(||) extremize the 



14 



free energy, 



V 



2v^7re \ Nf^ 



3 ^ 3v^ J 
Minimizing the free energy with respect to 5*, we obtain 

2/3 



(53) 



= (127rS)V3 I 1 + 



Nf 



(54) 



with the surface area, 



2/3 



1 + 



V 



2/3 



(55) 



Thus, we cannot tell, within the leading order of the thin-wall expansion, whether the F-ball 
is stable against the classical deformation and fragmentation into smaller ones. 

We then keep the terms up to the next-to-leading order in the thin-wall expansion to 
have 

v; = v;(°) + v;(i) , (56) 

with V}^^ smaller than V}°^ by the order of 5sc/R [H]. From Eqs.(||) and (0), K^^^ satisfies 



with the boundary condition. 



(0) 



V}^^ for \w\ > Ss 

The solution to Eq.(0) is given by 

ya) ^ __V6n\w\fl , l^cosh(^) 



(57) 



(58) 



2ew \Ri R2J smh^{^^) 



Iwl -|- 5s 



A 



-/ 



T 



where 



1 2 1 

f(x) = - sinh X cosh x H — tanh x sinh x — x 

•'^ ^ 3 3 3 



(59) 



(60) 



(see Fig.|^ for the illustration of V^^^ and Vj^^). Substituting the solutions, Vj'^^ and V}^\ 
into Fsc and keeping the terms up to the next-to- leading order, we get, in a manner similar 
to the derivation of Eq.(pl|), 

Fsc = + , (61) 
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W / JlT 



FIG. 5: The perturbative solutions for the electron potential energy VJ=. Here, Vq is the leading term 
and Vi the next-to-leading term in the thin- wall expansion. Note that Vq is rescaled as Vo/2ttT 
with T temperature, and Vi as eVi/\/^Tr with Vi = 2RiR2Vi/{Ri + R2) (Ri and R2 are the radii 
of principal curvature of the F-ball surface). 



where 29 



J |iXiri2| J \iXi -ft2 / 1X11X2 

Here C^c and Dgc are given by 

~ -50A^/% , (63) 

~ 25A^/% , (64) 

with Ve^^ = 2RiR2Vj^^ /{Ri + -R2)- Note that Dgc is positive in contrast to the case with 
the neutral F-ball where Dn in Eq.(p^ is negative. This feature is very important for the 
discussion on the stability below. From Eqs.(^) and (p^), we obtain 
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Let us take the curvature radius of the F-ball surface positive, and consider the fragmen- 
tation of the spherical F-ball into smaller F-balls. The negative sign of the first term in 
Eq.(|65D shows that the fragmentation is energetically promoted, which is the same as the 
neutral F-ball in Sec.|I|. On the other hand, the sign of the second term is positive for 

A > 1.2 X 10-^ , (66) 

and there should be the energy barrier to suppress the fragmentation of the F-ball. In this 
sense, electrically charged F-balls which are screened by electrons are metastable at 
T < O.lf (see the discussion in Appendix 0). 

The lifetime of the screened F-ball is estimated in Appendix |^ and is found extremely 
long, which makes the F-ball effectively stable. Thus, the screened charged F-balls can be 
the candidates for CDM. In order for them to be a main component of CDM, they should 
satisfy the constraints to be discussed in the next section. 

V. CONSTRAINTS ON THE PARAMETERS OF F-BALL 

We consider the phase transition of the universe where the biased Z2 symmetry is spon- 
taneously broken and the network of domain walls is formed [0. The number density of 



domain walls is soon diluted by annihilation to become ~ l/i?^, where H is the Hubble 
constant. Using the energy density of the domain walls, pyjoii ~ Sif, we get the ratio of 
Pwaii to the total energy density in the universe ptot-, 

Ptot ~ g.T^ ~ HMl, ' ^^'^ 
with (7* the effective degree of freedom |l^ and Mpi the Planck mass. The ratio will be of 
the order of unity when H ~ H^q = T^/M^i. 

The small violation of the Z2 symmetry gives an energy of the false vacuum bubble 
proportional to its volume. This volume energy increases with the universe expansion and 
gets comparable to the surface energy at ~ Hf = e/S. The domain walls then stop to 
expand getting round to produce F-balls. From the cosmological viewpoint, F-balls should 
be produced before the domain walls dominate the total energy density in the universe, 
since otherwise the energy of the universe would be dominated by blackholes made with the 
domain walls [12|. This condition is Hf > H^g, which is rewritten in terms of e. 



e > OAXv^M-i^ . (68) 



'pi 
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Let us consider the case where the F-balls are metastable and may contribute sizably to 
CDM in the present universe. As mentioned in Sec.^ we are not interested in the neutral 
F-ball, which is very tiny after the fragmentations, since we need too many new assumptions 
in introducing such an object and have only indirect means to detect it. We thus consider 
the charged F-balls. Since they interact with ordinary matter through the electric force, they 
would be easily observed in various experiments for the dark matter search if there exists 
enough cosmic flux. If the F-ball has a considerable contribution to CDM, the number 
density and its cosmic flux should be inversely proportional to its mass, 1/Mf. Since such 
events have not been observed so far, we should have Mf > 10^^ GeV p. This constraint is 
rewritten in terms of iVj, 



If the size of the F-ball is smaller than the typical screening length, its electric charge is not 
screened at the close neighborhood of the F-ball but screened far outside of it. In this case. 



the size of the unscreened F-ball R ~ 0.5yWj^(l + A/M^ff'^/X^/^ and the screenmg 

length 5sc ~ which gives a very small number of the fermion on the unscreened 

F-ball, Nf < 1000, from R < 5sc- This is incompatible with Eq.(]69D, and we see that the 
unscreened F-ball cannot give a sizable contribution to CDM. 

If the size of the F-ball is much larger than the screening length (i.e., R ^ (5sc), its electric 
charge is strongly screened as was discussed in the previous section. In this case, Eqs. (pT]) 
and (^) give Nj ^ 7000, which is compatible with Eq. (|69D . This allows screened charged 
F-balls to be candidates for the main component of CDM in the present universe. 

In the previous section, we found that the second term in Eq.(^) is crucial to enhanc- 
ing the stability and suppressing the fragmentation of the screened charged F-ball for the 
parameter satisfying Eq . (|66D . If we add the volume energy Ey ~ eV which has the desta- 
bilizing effect [0, we still have the metastable F-ball as far as the second term in Eq.(|65D 
overwhelms E^. Using Eq.(0), we express this condition as 




(69) 



the discussion on the stability made in Sec.|T| is valid. From Eqs. (|4^) and (0), we have 





(70) 
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VI. SUMMARY AND DISCUSSION 

We have considered the F-balls produced in the early universe due to the spontaneous 
breakdown of the biased Z2 symmetry. We have investigated the stabihty of the electri- 
cally neutral F-balls and the charged F-balls in order to examine whether they can sizably 
contribute to CDM or not. 

In the case of neutral F-balls, we have taken into account the correction to the thin-wall 
approximation up to the next-to-leading order, and found this correction plays an important 
role to enhance deformation and fragmentation of the F-ball into tiny thick-wall F-balls. If 
these tiny F-balls are a main component of CDM, the dark matter search experiments allow 
the F-ball to have only a weak interaction with the ordinary matter. We are not interested 
in such a new kind of weakly interacting particles since we need too many new assumptions 
in introducing new tiny particles and have only indirect means to detect such particles 

In the case of charged F-balls, two cases are considered: 

case 1) the size of a charged F-ball being smaller than the typical screening length. Such 
F-balls are found to have rather light masses, and they have a large number density in 
the universe if they are a main component of CDM. However, since the cosmic flux and 
the number density of the charged F-balls are observed to be very small, it is difficult for 
them to be a main component of CDM. (In this case, the electric charge of the F-balls is 
effectively unscreened. We have found that they decay into smaller F-balls through the 
tunneling effect though they are classically stable. Since this decay further increases the 
large number density of the F-balls, it merely decreases the possibility for them to be a 
main component of CDM.) 

case 2) the size of a charged F-ball being larger than the typical screening length. In this 
case, the electric charge of the F-ball is screened in the vicinity of its surface. Even though 
the long-range character of the Coulomb force is lost, the F-balls can still be classically stable 
in certain region of the parameters. This is due to the curvature effect that is obtained from 
the next-to-leading order correction to the thin-wall approximation. Such F-balls can be 
candidates for a main part of CDM in the present universe. In such a case, we have obtained 
constraints on the physical parameters: A > 1.2 x 10"^ ^AXv^M:^^ < e < 100\'^^%*Nf ' 
and f < 3 X 10^ GeV. It is interesting to note that the symmetry breaking scale v constrained 
above is not too far from Electroweak or SUSY breaking scale. These above constraints help 
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us to make a realistic model of the F-ball. 

Recently the dark matter distributions of galaxies and clusters were observed and com- 
pared with the dark matter simulations. Some find consistency of the collisionless CDM 
simulation with the observations |T^, |T3], while others find inconsistency and propose the 
self-interacting dark matter (SIDM) with the ratio of the cross section to the mass as large 
as 10~^^ — lO'^'^cm^ GeV"^ In the case of screened charged F-ball, the geometrical cross 
section per unit mass is not so large, 

^ ^ X 10-28 cm^GeV-i , (72) 

and, moreover, the typical momentum transfer squared ~ ~ "^^/^z is too small to give 
a large angle scattering which appears to be necessary for smoothing the Halo central density 
profile. Thus, the present model of the F-ball is not adequate for SIDM. We, however, do not 
discuss this issue further in the present paper since the above problem is still controversial 

EMM- 



All through the present paper, we have dealt with the heavy fermions as being distributed 
two-dimensionally on the F-ball surface, and neglected the spreading of the fermion in the 
direction normal to the surface. This spreading may have some contributions to the curvature 
effect energy of the F-ball. We will discuss this issue elsewhere. 

APPENDIX A: VOLUME ELEMENT IN NEW COORDINATES 

We introduce three local orthogonal coordinates {u, v, w) for the F-ball, the origin of 
which is on the surface of an F-ball, and derive a representation for the volume element in 
this frame. The coordinates, u and v, represent the F-ball surface, x = x/(M,f), and are 
taken along the lines of curvature |^ on the surface. We define the first and the second 



fundamental forms, 

/ = E{u, v)du^ + 2F{u, v)dudv + G{u, v)dv^ , (Al) 
J = L{u,v)du^ + 2M{u,v)dudv + N{u,v)dv^ , (A2) 

with the coefficients, 

E{U,V) = Xf,u^f,u , L{U,V) = Xf^uu^ , 
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F{u,v) = x/,„x/,„ , M{u,v) = x/,„„n , 

G(u,t;) = x/^^x/_„ , N{u,v) = :s.f^yyn , (A3) 

where the subscripts of u and v stand for the partial derivative with respect to them, and 
n{u, v) is a unit vector normal to the surface defined as 

niu.v) = -r^ . 

Since we take u =const. and v =const. as the lines of curvature, we obtain 



(A4) 



F{u,v) = M{u,v) = . (A5) 
We define the third coordinate w by 



:k{u,v,w) = :Kf{u,v) + wili{u,v) , (A6) 

taking w positive (negative) for the outside (inside) of the F-ball. We take w = at the 
place in the wall where the scalar field vanishes. 
Using the orthogonal relations, 

x/,„x/,„ = x/,„n = x/,^n = ndn = , (A7) 

we express the line element squared, 

d^x = {-Kf^udu + Xf^ydv + ndw + wdn)"^ 

= Edu^ + Gdv^ + 2w(x/^^,dndM + Xf-^dndv) + w'^dn^ + dw^ . (A8) 

From Eg. (|A7|) and M{u,v) = in Eg. (|A5|) , we obtain 

n« = -^^/,« ' (A9) 
N , 



Thus, we obtain 



xjudn = xj„ (n^d-u + nt,df ) = —Ldu , (All) 
x/^^dn = x/_„ (n„dM + w^dv) = —Ndv , (A12) 



dn^ = (liudu + n^dv) 

= + G 
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du^ + —dv^ . (A13) 



Substituting Eqs. (|XTTl) to (|M3|) into Eq.(^, we get 



d'x = E 1 



wL 



du' + G[l- 



wN 



E J V G 

Calculating a determinant of the line element squared, we obtain the volume element 



(A14) 



d'^x 



'EG 1- 



wL 



1 - 



wN 



E J \ G 

wN 



dudvdw 



dSdw 



where we use the relation, 



d^ = v^dudv . 



(A15) 



(A16) 



Taking Ri and R2 as the radii of principal curvature of the surface, we obtain the Gaussian 
curvature, 

K 



1 



LN 



R1R2 EG 



the mean curvature. 



and the volume element. 



EN + GL 
2EG 



d^x 



-Kl -K2 



dS'diy . 



(A17) 
(A18) 

(A19) 



APPENDIX B: CURVATURE EFFECT IN NEUTRAL CASE 



We here derive Eq. (|2T|) to (p^) in Sec.|T|, assuming Ri, R2 > 0. Noting that higher 
order effects in Etot arise only from Eg, we estimate Eg keeping the terms up to the order of 
{5n/RY- Substituting Eq.(|18D into Eq.(^, we obtain 



where 



Eg = + E« + ) 



dS 



dw 



2\2 



W W \ i 1 ; 2 A , ,2 .)\ 



A 



+ \ 0001 + 20o(0o - ^^)0i 



(Bl) 



(B2) 



(B3) 
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) = fds dw 



w"^ f 1 ; 2 A 



R1R2 12 8 



+-(0o-^Y 



Substituting Eq.(|T5D into Eq.Q, we get 



=T. Us 



with S = 2-\/Ai^'^/3. The integration of Eg. ( |B3| ) vanishes, 



E^}^ = 



because the integrand is an odd function of w. Using the relation, 0o 
derived from Eg. ([T3|) , we express the first hne in the r.h.s. of Eg. (P^) as 

I- /■+CXD _ 

(1st hne) = dS dw 



-R1-R2 



The second hne is rewritten as 



(2nd hne) = dS dw 



+00 



w ^ w 
Ri R2 
1 



-00 + ^0o(0O - 



Ri R2' 
+00 /I X 

dS I dw [ — + — ] 0ov^i 
00 \ rii ri2 . 



From Eg.(^), we get for the third hne. 



(3rd hne) = j d5 
1 



^ f 1 •■ A 1 



|-+oo /I 1 

z J J —00 \JXi 1x2, 



From Eg.(|B3) to (§|), we obtain 

J J —00 I 111 

Using the dimensionless guantity. 



w^ ■ 2 1/1 1 



1x2 ^ \ -Ki 1x2 , 



2R1R2 

Ri + R2 
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we express E^^^ as 



dS 



+00 



dw 



w 



1 / 1 



R1R2 4 \Ri 

CnfdS-^ + oJdS (-1-^ 
J 1X11X2 J \1X\ -it; 



1 
'R2 



with 



Cn 



4 j-00 

From Eq.(0) with Eq.(|2g), we get 

^^sech^ ( -r-\ ^cosh^ ( — 



+00 / .2 . ~ 

dw wVo -0001 

—00 \ 

1 f+°° ■ ~ 

dw 0001 . 



(B12) 

(B13) 
(B14) 



h[w) 



— ^smh — I tanh ( — 



H — ^^sech^ 



, /4w\ , /2w\ 12w 

smh — + 8 smh — H — - — 



with 6n = 2/y/Xv. Substituting Eqs.(||) and (|BT5| ) into Eqs.( pT3|) and ( pl^) , we 

C'n, 



(B15) 
obtain 



00 r 2 

dw |4iy^sech^ii; — 2sech^w + 2 — -tanh'^w 
— -sech'^w(sinh4iu + 8sinh2w + 12iD)| 



/ dw Aw sech w ~ 1=^ 



(B16) 



5f 



dw 



■-sech^w H tanh^w 

2 2 6 



— sech^w(sinh4w + 8sinh2M} + 12w7)| 



18v^ ■ 



(B17) 



APPENDIX C: ENERGY OF UNSCREENED F-BALL 



We derive the thin-wall representations for the energy of an unscreened charged F-ball, 
Eqs.([3T|), (|3^ ) and (|33D . We first consider the surface energy, 



Es = j: I dS . 



Using Eq.(^), we get the metric of the surface, 

dx^ = (r^ + f2) de^ + r^m'^edip'^ 



(CI) 



(C2) 



24 



with the symbol 'dot' standing for the derivative with respect to 9. From Eq.i 



d^ = r^yi+y d^]. 

Substituting dS in Eg. ( |C3|) into Eq.( |Cl|) , with r expressed in Eg. (p8[) , and usin^ 
orthogonahty relations, 

1 



4-k' 



2/ + 1 

2/ + 1 



we obtain Eq.(|3l|) with Eq.(|3|). 

We next consider the fermi energy, 



Using Eg. (pO|) and Eg. (P3|) , we rewrite Eg.(P6D as 



3/2 



3R J Atc (1 + 55)1/2 
2nT rdQ 



3R 



where 

^daj + -60, 

This gives Eg.(|32|) with Eg.®. 

We finally consider the Coulomb energy 

.2 



6F = -6a f + -6a^f — 6r + 5r^ — -5r^ — SatSr . 



"I"! 



X — X 



Using the expression, 



47r|x — x' 



d'^k 



we rewrite Eg.(p9D as 

47r^ Jq J An 
Using the following expansion relation 



' dk 



(27r)3A;2 
dQx 



^ik(x-x') 



Air 



{l + 6af)e 



ikx 



kx ) 



1=0 

oo 



1=0 m 



2|), we obtain 
(C3) 

the following 

(C4) 
(C5) 



(C6) 



(C7) 



(C8) 



(C9) 



(CIO) 



(Cll) 



(C12) 
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(here ji{kr) and Yim{^) are the spherical Bessel function and the spherical harmonics, re- 
spectively, and Qkx is the angle between k and x), and the orthogonality relation, 

J dfifcY'im(fifc)yj,*„,(fifc) = 5ii'5mm' , (C13) 

we obtain 

E^ = ^J^y-l^ , (CM) 

47r2 ^ 2/ + 1 ' ^ ^ 

where 

xji{kr)3i{kr')Pi{cose)Pi{cose') . (C15) 

Here, we have introduced the procedure, lim^^O) in order to safely expand the spherical 
Bessel function in the integrand, 

jiikr) ~ jiikR) + 5rkR3[{kR) + Ur^kRfj'^ikR) . (C16) 

Keeping the terms up to the second order, we obtain 

D« = {2l + lf\\mj^dke-'^j'^^^^Pi{cosd) 

.J^^P^{cose')l['\e,9'), (C17) 



where 



= {l + 6(Tf + 6(r'f + 6af6a'f)ji{kRf 

+ (1 + 5(7/ + 6a'f){6r + 5r')kR]i{kR)][{kR) 
+5r5r\kRf{j[{kR)f + 6r5r' {kRf{j[{kR)f 

+ 1(5^2 ^ 5/2)(^^)2^.^(^^)_^..(^^) _ (C18) 

Integrating Eg. ( |C17] ) with respect to the angular variables, we get 

= lim r dke-"^I^2^ , (C19) 

where 

I? = {^m + c^jK^^)' + (25^0 E + 2a,Q^ kRji{kR)Jl{kR) 



n= 



^ al 



+4{kRf{j[{kR)f + 5,0 E 7^(^i?)'jo(^i?)jo'(A:i?) • 

^1 2n + 1 



(C20) 
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Using the following integration formulae, 

f d.e-0,(*fi))^ = ^^^ + O(.), (C21) 

dke-^^kR3i{kRmkR) = -^^^—^ + 0{e) , (C22) 







^ dke-^\kRnjlikR)r ^ '"^^^ + ^^ 



2e 2(2/ + l)i? 

+0(e, or e loge) , (C23) 

dke-^\kRfUkR)f^{kR) = -- + — 

+0(e, or e loge) , (C24) 

we obtain Eq.® with Eq.®. 

APPENDIX D: DECAY RATE OF UNSCREENED F-BALL 

Unscreened charged F-balls are metastable and fragment through the quantum tunneling 
(see Sec.|T| for the details). Here, we roughly estimate the rate of fragmentation, assuming 
a very simplified model of fragmentation after the ellipsoidal deformation from the spherical 
shape 10 . 



We express the surface of the ellipsoidal F-ball as 

x'^ iP' 

with 7] being a constant which we call 'deformation parameter'. Note that R depends on 77 
since it is determined as to minimize the energy of the F-ball. The surface is also expressed 
in the polar coordinates as 

r = Rf{e) = . (D2) 

d cos^d + ?72sin 6 

The scalar field 0(r, 6) should satisfy the Euler-Lagrange equation to be derived from Eq.(|), 

V^</> = 0" + h' + -,^+ = ^4>{4>' - v') , (D3) 

r r^tanfc' 2 

with the symbols, 'dash' and 'dot', standing for d/dr and d/d9, respectively. Here, the 
boundary conditions are imposed, 

p{r 00) = +v 

P{r = Rf) = (D4) 
(j){r — s> 0) ~ —V . 
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(The last one becomes exact when Rf tends to infinity (see Eq.( |D7|) ).) Using the following 
variable, 

f = r-Rf{e), (D5) 



and keeping the terms of the leading order for Rf oo, we can rewrite Eq. (|D3|) as 



1 ^./.2 2 



02 Qf2 2 



- v') , (D6) 



to give the solution, 



with 5n = 2/ \/\v and 



i;tanh|G^^^y^| , (D7) 



1 /8/?/N2 
32 I ^2 a2 



cos + 7]^ sin 6^^ 



(D9) 



Vcos6'2 + r]^ sin 6*2 

We regard the deformation parameter rj as the time-dependent dynamical variable. Sub- 
stituting into the Lagrangian density of Eq.(^ with the Fermi, the Coulomb and the 
volume effect added to it, and integrating it over the whole space, we get the Lagrangian, 

L = hi{r^)rf - W{v) , (DIO) 
where /i(?7) and W{r]) are estimated as follows. First, h{ri) is expressed as 

where 

N, = f dx x^^'^-i) {(r/2 - l)x^ - ry2}-'^ {^4 _ (^4 _ i)a;2|-i/2 ^ 2, 3) . (D13) 



Since 1/cosh {6(r — in the integrand in Eq. (pil|) is negligibly small except at 



around r ~ Rf, we approximate the quantities which are smooth at r ~ i?/ to those at 
r = Rf. Now let us estimate W{ri) = + Eg + Ef + Ec. We easily obtain the volume 
energy, 

= -^e . (D14) 
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In order to estimate Es, Ef and E^., we introduce new variables 9 and (p defined by 



X = R sin 6 cos (p 
y = R sin 9 sin (p 
z = TjR COS 9 . 



(D15) 



The surface element is expressed as 



d^ = R'^^cosW + r]^sm'^9sm9d9d^ 
= R^ J k{u)dud(p , 



(D16) 



where u = cos 9 and k{u) = u"^ + ri^{l — v^). From Eg. ( pi 61) , we obtain the surface energy, 

Es = T. j dS = A-kHR^ du ^fk{u) . (D17) 



We express the fermion number density, 



l+Y.niPi{u) 



(D18) 



AnR^^k{u) \ 1=1 / 

where Pi{u) is the /'th order Legendre polynomial. It is easy to see that the above expres- 
sion satisfies the condition for the fermion number Nf to be conserved (see Eq.(^). From 
Eg. ( pi 61) and Eq.( pi8| ), we obtain the fermi energy. 



3 J ^ 3R Jo k{u 
Let us estimate the Coulomb energy, 

K='-l IdSdS'.^. 



1 1/4 



57r 



X — x 



Using the relation. 



47r|x — x'l 



^iq(x-x') 



(27r)3g2 



d^q 



{2Tif(tk{u^) 



3*q(x-x') 



with 



X = (x^, Xy, -x^) = -R(sin^cos sin^sin cos ^) , 



q = {(lx,qy,mz) 



and Ug = Qz/q, we get 



.2 Ar2 



E. 



Nfri 



47r^ Jo 



dq 



Ank{uq) 



dQ 



Atx 



\ l + Y.niP,{u)]e 



iqx 



1=1 



(D19) 

(D20) 

(D21) 

(D22) 
(D23) 

(D24) 
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With the aid of 

3i{(iR)ji'{qR) 



TT 



for I — I' = even integer 



(D25) 



10 2(21 + 1)R' 

the integration in Eq. ( |D24|) can be carried out in the same way as in Appendix to give 

^1 



S-kR Jo k{u) 



E 

1=1 



21 + 1 



(D26) 



We then minimize W{r]) with respect to R and rii. As expected from the perturbative 



analysis in Sec. Ill, we verify numerically that W{ri) takes the local minimum value Wmin 
at ?7 = 1 in a certain range of the parameters (see Fig.^ for the case, e = 0, A = 1, and 
Nf = 1000). In the following, we consider such an F-ball which has the local minimum 




Tl 



FIG. 6: The potential is displayed as a function of rj in case the F-ball is not screened for the 
parameters shown in the graph. There is an energy barrier, which is comparable to v. The 
distribution of the fermion is shown in Fig.^. 

value of W at 1] = 1, that is, is stable against the deformation from the spherical shape. 
We finally proceed to estimate the decay rate of the metastable F-ball, using the Euclidean 



action method in Ref.|jT^. We express the Euchdean action as 

A = Jdt [Ihir^W + Viv)] , (D27) 
where V{r]) = VF(?7) — Wmin ("we note V{1) = 0). Defining the following variable, 

^ = I'dr^'^/hi^ , (D28) 
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FIG. 7: The number density of the heavy fermion as a function of 6 for r] = 0.8, 1, 1.7. The density 
is constant for ry = 1. It gets large at the equator for rj = 0.8 and at poles for rj = 1.7. They are 
graphically displayed in Fig.^. 

we express the action as 

A=jdt [le+no] , (D29) 

where the potential V{^) = V{ri{^)) has a local minimum = at = 0. According to 



Ref . |117| , the decay rate F is estimated semi-classically, 

T = n\K\ e--^o/^ , (D30) 

where is the classical action of the bounce solution and K is the prefactor, both of which 
are evaluated below. The bounce solution for ^ should satisfy the following equation, 

^ = -e + V'iO = , (D31) 

with the boundary condition, 

^ for t ^ ±oo . (D32) 

From Eq.( P31|) , we get = 2V{^) for the classical solution, and obtain the classical Eu- 



chdean action. 



= 2£dv ^/2h{r])V{v) , (D33) 
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(a) 



(b) (c) 



FIG. 8: The distribution of the heavy fermions in case the F-ball is not screened, (a) The F-ball 
becomes cigar-Hke shape for ij > 1 and the fermions gather at poles due to the repulsive electric 
force, (b) The F-ball is a sphere for ij = 1 and the distribution is constant. (c)The F-ball has the 
pancake- like shape for ij < 1, and the fermions gather at the equator due to the repulsive force. 

where rjc {rjc > 1) and (^c > 0) are defined as V{ric) = V{^c) = 0. The constant K is 
defined as 



rj7 det{-£, + v (oy 



where the symbol 'det" means the determinant in which the zero-mode contribution is 
removed [1^]. In our case, K is calculated as 



3/4, , m~j2v{o 



h [ Ml) ) 



Taking parameters v = 10^ GeV, Nf = 1000, and A = 1 for example, we obtain — 
550. Thus, the lifetime of the unscreened charged F-ball is much longer than the age of 
the universe in this case. (If the production temperature of the F-ball Tf is as high as 
Tf ~ Tph ~ V, the F-ball would decay through the thermal fluctuation. We here consider 
the temperature Tf ^ O.lv Tph where the lifetime of the F-ball is found to be long enough 
to let it survive till present.) 
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APPENDIX E: DECAY RATE OF SCREENED F-BALL 



Here, we roughly estimate the decay rate of the screened F-ball, assuming the very 
simphfied model, in which the surface has the following shape, 

Rf = Ro{l + a2P2{cose)), (El) 

where P2 is the 2nd order Legendre polynomial and Rq and 02 are arbitrary constants. This 
allows the F-ball to have a concave shape. In the following we consider the F-ball energy 
as the function of the area S and the parameter 02. Then, the parameter Rq is a function 
of S and 02. We use the same technique as that in Appendix 0. We here regard 02 as the 
time-dependent dynamical variable instead of 77 in Appendix 0. All we need to obtain are 
h{a2) and V{a2) in the Lagrangian, 

L = h{a2)d2' - V{a2) . (E2) 

We can easily obtain h{a2) by replacing Rf in Eqs. (|D8| ) and (p 1 1|) with that defined by 
Eq.Q. 

The potential V{a2) can be obtained from Eq . (|65D since fI;-^^ depends only on 5* and 
not on the shape of the F-ball. We calculate the following quantities using the method 
introduced in Appendix ^ 



Xc ^ dS —— , (E3) 

J |/iiri2| 

Xd = [ dS (- 1 y \RiR2\ ^^^^ 

J \Ri R2J R1R2 

The surface is expressed as 

X = x/-(6', = -R/(^^)(sin6' cos V9, sin6' sin y9, cos6') , (E5) 
in the polar coordinates. The first and the second forms are 

/ = E{e, 0)d^2 + 2F{e, (j))ded(j) + G{e, 0)d02 , (eg) 

J = L{e, (f))de^ + 2M{e, 0)dM0 + a^(^, 0)d0^ , (et) 



with 



m <P) = R}+ Rj, m 0) = -j=^^{ii'/R} - 2i?? - R}), 

F{9,4>) = 0, M(^,0)=O, 

G{9, 0) = R} sin^ 9 , N{9, 0) = sin 9{-Rf sin 9 + Rjcos 9) , (E8) 

V ^/ + ^/ 

33 



where the symbol 'dash' denotes the derivative with respect to 6. From these, the Gaussian 
curvature in Eq. ([A17|) and the mean curvature in Eq. (|A18D are expressed in terms oi Rf. 
We then obtain 



1 



R1R2 



K 



{R!f sin e - Rf cos e){-R'f'^R} + 2RJ + R 



Rf smd{R) + R:ff 



(E9) 



Rl i?2 



{EN - GLf 

[{R) + R'f){-R'f sin e + Rf cos 6) + {-R'j'^R} + 2RJ + Rj) Rf sin O}' 



R)sin^e{R) + R:ff 



(ElO) 



From these expressions, Xc and Xo are written as 



Xc = 47r / d^ 



(i?^ sin 9 - Rf cos e){-R'f^R} + 2R^f + R 



(Ell) 



Xd = An I de 



{Rj + Rjr^ 

[{Rj + R'f){-R'f sine + Rf cos 6) + {-R'f^Rj + 2RJ + R))Rf smO'f \k 



Rf sin e{R} + Rf 



K 

(E12) 



We note that these quantities do not exphcitly depend on the parameter Rq or S and only 
depend on a2. Using the free energy. 



(E13) 



we evaluate the potential. 



7(2), 

tot ' 



02) 



(E14) 



We show the potential as a function of a2 in Fig.^ for the case T = O.lf. The cusp in 
Fig.|| arises because of the non-analytic feature of Xc and Xd. At the top of the cusp, the 
Gaussian curvature becomes zero at equator and turns to be negative for larger 02. This 
is illustrated in Fig.^. One can observe that the curvature at the origin is positive and 
the potential becomes negative for 02 larger than a2c- These facts mean that the spherical 



F-ball is metastable as discussed in Sec.|I^. Because the energy barrier is much larger than 
the temperature, the decay of the F-ball due to thermal fluctuation can be neglected. We 
next consider the decay due to the tunneling effect. 
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FIG. 9: The potential of the screened F-ball as a function of a2. The curvature at the origin 
is positive and the potential becomes negative for a2 larger than a2c at T = 0.1^;. These facts 
mean that the spherical F-ball (02 = 0) is metastable. The energy barrier is much larger than the 
temperature. This tells us that the decay due to the thermal fluctuation can be neglected. 

We consider the classical Euclidean action, 

Ao = 2 das j2h{a2)V{a2) . (E15) 
Jo 

We can apply Eq.( |Dll| ) to obtain h{a2) and see its magnitude is 0(S/2'^). We find that the 
barrier height is O(100f ) and a2c is 0(1) from Fig.^. From these, the classical Euclidean 
action is roughly estimated as Aq ^ Nf. Since the number of the heavy fermion is large, 
Nf > 10^^, from Eqs.(^) and (O); ^^e classical Euclidean action is extremely large and the 
screened F-ball is effectively stable. Although we use the very simple deformation model 
and take the special case of T = O.lf , we think that the result of very long lifetime is rather 
general for the temperature smaller than v since the result is only due to the macroscopic 
property of the F-ball. 
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becomes slightly non-convex at the equator for 02 = a2c ~ 0.32 and more narrow in the middle for 
a2 = 0.4. 
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